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ABSTRACT 

We consider the recently estimated corrections Au> to the Newto- 
nian/Einsteinian secular precessions of the longitudes of perihelia vj of several 
planets of the Solar System in order to evaluate whether they are compatible with 
the predicted precessions due to models of long-range modified gravity put forth 
to account for certain features of the rotation curves of galaxies without resorting 
to dark matter. In particular, we consider a logarithmic-type correction and a 
f(R) inspired power- law modification of the Newtonian gravitational potential. 
The results obtained by taking the ratio of the apsidal rates for different pairs 
of planets show that the modifications of the Newtonian potentials examined in 
this paper are not compatible with the secular extra-precessions of the perihelia 
of the Solar System's planets estimated by E.V. Pitjeva as solve-for parameters 
processing almost one century of data with the latest EPM ephemerides. 

Subject headings: Experimental tests of gravitational theories; Modified theories of 
gravity; Celestial mechanics; Orbit determination and improvement; Ephemerides, 
almanacs, and calendars 
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Introduction 



Dark matter and dark energy are, possibly, the most severe theoretical issues that 
modern astrophysics and cosmology have to face, since the available observations seem 
to question the model of gravitational interaction on a scale larger than the Solar 
System. In fact, the data c oming from the galactic rotation curves of spiral galaxies 



( jBinney and Tremaine 



19871 ) cannot be explained on the basis of Newtonian gravity or 
General Relativity (GR) if one does not introduce invisible dark matter compensating for 
the obs e rved i nconsistency of the Newtonian dynamics of stars. Likewise, since 1933, when 



Zwickyl (119331 ) studied the velocity dispersion in the Coma cluster, there is the common 



agreement on the fact that the dyn amics o 



basis of Newtonian gravity or GR ([Peebles 



cluster of galaxies i s poorly understood on the 



1993 



Peacock 



19991 ) if the dark matter is not 



taken into account. In order to reconcile theoretical models with observations, the existence 
of a peculiar form of matter is postulated, the so called dark matter, which is supposed to 
be a cold and pressureless medium, whose distribution is that of a spherical halo around the 
galaxies. On the other hand, a lot of observations, such as the light cur yes of the type la 



supernovae and the cosmic microwave background (CMB) e xperiments (IRiess et al 



1998 



Perlmutter et al. 



1999 



Tonry et al. 



2003 



Bennet et al. 



20031 ) , firmly state that our Universe 



is now undergoing a phase of accelerated expansion. Actually, the present acceleration of 
the Universe cannot be explained, within GR, unless the existence of a cosmic fluid having 
exotic equation of state is postulated: the so-called dark energy. 

The main problem dark matter and dark energy bring about is understanding their 
nature, since they are introduced as ad hoc gravity sources in a well defined model of 
gravity, i.e. GR (or Newtonian gravity). Of course, another possibility exists: GR (and 
its approximation, Newtonian gravity) is unfit to deal with gravitational interaction at 
galactic, intergalactic and cosmo logical scales. The latter viewpoint, led to the introduction 
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of various modifications of gravity (INojiri and Odintsovl 120071 ; lEsposito-Faresd 120081 ) 



However, we must remember that, even though there are problems with galaxies and 
clusters dynamics and cosmolo gical observ atio ns, GR is in exc ellent agreement with the 



Solar System experiments (see IWilll (j2006l ) and iDamourl (j2006l )): hence, every theory that 
aims at explaining the large scale dynamics and the accelerated expansion of the Universe, 
should reproduce GR at the Solar System scale, i.e. in a suitable weak field limit. So, the 
viability of these modified theories of gravity at Solar System scale should be studied with 
great care. 

In this paper we wish to quantitatively deal with such a problem by means of the 



secular (i.e. averaged over one orbital revolution) precess ions of the 



perihelia w of some planets of the Solar System (see also (jlorio 



2007a 



ongitudes of the 



bj)) in the following 



way. Generally speaking, let LRMOG (Long-Range Modified Model of Gravity) be a given 
exotic model of modified gravity parameterized in terms of, say, K, in a such a way that 
K = would imply no modifications of gravity at all. Let P(LRMOG) be the prediction of 
a certain effect induced by such a model like, e.g., the secular precession of the perihelion 
of a planet. For all the exotic models considered it turns out that 



P(LRMOG) = Kg(a,e), 



(1) 



where g is a function of the system's orbital parameters a (semimajor axis) and e 
(eccentricity); such g is a peculiar consequence of the model LRMOG (and of all 
other models of its class with the same spatial variability). Now, let us take the ratio 
of "P(LRMOG) for two different systems A and B, e.g. two Solar System's planets: 
7 ? a(LRMOG)/'Pb (LRMOG) = #a/<?b- The model's parameter K is now disappeared, 
but we still have a prediction that retains a peculiar signature of that model, i.e. g^/ 9b- 
Of course, such a prediction is valid if we assume K is not zero, which is just the case 
both theoretically (LRMOG is such that should K be zero, no modifications of gravity at 
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all occurred) and observationally because K is usually determined by other independent 
long-range astrophysical/cosmological observations. Otherwise, one would have the 
meaningless prediction 0/0. The case K = (or K < K) can be, in stead , usually tested by 



taking one perihelion precession at a time, as already done, e.g., by 



Ioriol (2007a 



If we have observational determinations O for A and B of the effect considered 
above such that they are affected alsqj by LRMOG (it is just the case for the purely 
phenomenologically estimated corrections to the standard Newton-Einstein perihelion 
precessions, since LRMOG has not been included in the dynamical force models of the 
ephemerides adjust ed to the plan et ary dat a in the least-square parameters' estimation 



process by Pitjeva (IPitjeva 



2005a 



2006bl )). we can construct Oa/O b and compare it with 
the prediction for it by LRMOG, i.e. with qa/qb- Note that 50/0 > 1 only means that 
O is compatible with zero, being possible a nonzero value smaller than SO. Thus, it is 
perfectly meaningful to construct Oa/O b . Its uncertainty will be conservatively evaluated 
as \1/0 B \80 A + \O a /O b \SO b . As a result, O a /O b will be compatible with zero. Now, the 



question is: Is it the same for qa) '<7b as well? If yes, i.e. if 

Oa _ P A (LRMOG) 
Ob~ ~ V B (LRMOG) 

within the errors, or, equivalently, if 

O a P A (LRMOG) 



(2) 







(3) 



O b V B (LRMOG) 

within the errors, LRMOG survives (and the use of the single perihelion precessions can be 
used to put upper bounds on K). Otherwise, LRMOG is ruled out. 



l \i they are differential quantities constructed by contrasting observations to predictions 
obtained by analytical force models of canonical effects, O are, in principle, affected also by 
the mismodelling in them. 
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The paper is organized as follows. In Section [2] we apply this approach to a test particle 
in motion around a central mass M whose Newtonian gravitational potential exhibits a 
logarithmic-type correction. Then, in Section [31 we consider the power-law modification of 
the gravitational potential inspired by f(R) extended theories of gravity. Comments and 
conclusions are outlined in Section HI 



2. The perihelion precession due to a 1/r force 



Logarithmic corrections to the Newtonian gravitational potential have 



jeen recently 



used to phenomenon 



1987 



2007 



Cretton 



Sobouti 



etal 



ogica l 



2000 



y tackle the problem of dark matter in ga l axies ( 



Kinnev and Brisudqva 



2003). E.g.. lFabris and Camposl (120071 ) used 



2001 



Kirillov 



2006 



van Moorsel 



Fabris and Campos 



VL = -aGMln 



r 



(4) 



where a has the dimension of L l , to satisfactorily fit the rotation curves of 10 spiral 
galaxies getting 

a « -0.1 kpc -1 . (5) 

The extra-potential of eq. (jlj) yields an additional 1/r radial force 

aGM . 



A 



in 



r. 



(6) 



Various theoretica. 
E.g., according to 



justifications have been found for a logarithmic-like extra-potential. 
Kirillov I2OO6 ). it would arise from la rge-scale discrep ancies of the 



topology of the actual Universe from the Fridman space; 



Sobouti 



(120071 ) obtained it in 



'or another ex ample of a 1/r extra-force and its connection with galaxy rotation curves 



see ( iSanders 



20061 ). 
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the framework of the f(R) modifications of general relativity getting preliminarily flat 



rotation curves, t 
of M OND, while 



re Tully-Fisher relat i on (a dmittedly with some reservations) and a versi on 



1995 



Fab ris and Campos 



Capozziello et al. 



(120071 ) pointed out that string-like objects (jSoleng 



20061 ) would yield a logarithmic-type potential with a related to the 



string tension. 

Alternative tests of such proposed correction to the Newtonian potential, independent 
of the dark matter issues themselves, would be, of course, highly desirable and could be, in 
principle, conducted in the Solar System. This is, indee d, considered one of the task s to be 
implemented in further investigations by 



Soboutil (120071 ): 



Fabris and Campos! (120071 ) argue 



that, in view of their extreme smallness due to eq. (151), no detectable effects induced by eq. 
(JHJ) would be possible at the level of Solar System. 



E.V. Pitjeva has recently processed almost one century of data of different types 
for the major bodies of the Solar System in the effort of continuo usly improving the 



EPM2004/EPM2006 planetary ephemerides (|Pitievall2005al . l2006bQ . Among other things, 
she also simultaneously e stimated corre ctions Am to the secular r ates of the lon gitudes of 



perihelia w of the inner (IPitjeva 



2005bl ) and of some of the outer (IPitjeva 



2006a 



bj) planets 



of the Solar System as fit-for parameters of a global solution in which she contrasted, in a 
least-square way, the observations to their predicted values computed with a complete set 
of dynamical force models including all the known Newtonian and Einsteinian features of 
motion. As a consequence, any un-modelled exotic force present in nature is, in principle, 
entirely accounted for by the obtained apsidal extra-precessionqfl Am. See Table [TJ for the 



3 Of course, since modelling is not perfect, in principle, Am include also the mismodelled 
parts of the standard Newtonian/Einsteinian effects, but it turns out that their inclusion in 
the following computation do not alter our conclusions. 
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inner planets and Table [2] for the gaseous giant ones. In regard to them, we must note that 
modern data sets cover at least one full orbital revolution only for Jupiter, Saturn and, barely, 
Uranus; this is why no secular extra-precessions of perihelia for Neptune and Pluto are today 
available. 

In order to make a direct comparison with them, we will now analytically work out the 
secular effects induced by the extra-acceleration of eq. © on the pericentre of a test particle. To 
this aim, we will treat eq. © as a small perturbation of the Newtonian monopole. The Gauss 
equation for the variation of m under the action of an entirely radial perturbing acceleration A r is 



dm 
~dt 



vT 



nae 



-A r cos /, 



(7) 



in which a is the semimajor axis, e is the eccentricity, n = yj GMf a 3 is the (unperturbed) 
Keplerian mean motion related to the orbital period P by n = 2n/P, and / is the true anomaly. 
After being evaluated onto the unperturbed Keplerian ellipse, eq. © must be inserted into eq. 
(J7J); then, the average over one orbital period P must be performed. It is useful to use the eccentric 
anomaly E by means of the relations 

r = a(l — e cos E), 



dt = ^- ec n osE) dE, 



(8) 



/( 
= I 



c os E—e 
-e cos E ' 



sin/ = 



e cos E 



On using 



it is possible to obtain 



2tt 



277 



cos E — e 
1 — e cos E 



dE 



(^)ln = -« 



GM{1 - e 2 ) / -1 + 



(9) 



(10) 
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' ?able 1: Inner planets. First row: estimated perihelion extra-precessions, from Table 3 of 



( IPitjeval l2005bl ). The quoted errors are not the formal ones but are realistic. The units 
are arcseconds per century (" cy _1 ). Secon d row: semima jor axes, in Astronomical Units 



(AU). Their formal errors are in Table IV of (jPitjevall2005al ). in m. Third row: eccentricities. 
Fourth row: orbital periods in years. 



Mercury 



Earth 



Mars 



(Aw) (" cy" 1 ) -0.0036 ± 0.0050 -0.0002 ± 0.0004 0.0001 ± 0.0005 

a (AU) 0.387 1.000 1.523 

e 0.2056 0.0167 0.0934 

P (yr) 0.24 1.00 1.88 



Table 2: Outer planets. First row: estimated perihelion extra-precessions ( lPitjevall2006al Jbl) . 
The quoted uncertainties are the formal, statistical errors re-scaled by a factor 10 in order to 
get the realistic ones. The units are arcseconds per century (" cy _1 ). Seco nd row: semim ajor 



axes, in Astronomical Units (AU). Their formal errors are in Table IV of (jPitjevall2005al ). in 
m. Third row: eccentricities. Fourth row: orbital periods in years. 



Jupiter 



Saturn 



Uranus 



(Aw) (" cy- 1 ) 0.0062 ±0.036 -0.92 ± 2.9 0.57 ± 13.0 

a (AU) 5.203 9.537 19.191 

e 0.0483 0.0541 0.0471 

P (yr) 11.86 29.45 84.07 
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Note that eq. (|10p is an exact result with respect to e. Eq. (|10p agrees with the precession 
ob tainable dividing; by P t he ad imensional perihelion shift per orbit A# p (log) worked out 



by. 



Adkins and McDonnell 



Adkins and McDonnell ( 



20071 ) within a different perturbative framework; note that for 



20071 ) the constant a has the dimensio ns o: ' 



-GMa must be performed in (50) of (jAdkins and McDonnell 



M L 2 T~ 2 , so that tha t 



2003) to 



the substitution a/m - 
retrieve eq. (fTUj) . 

It may be interesting to note that for the potential of eq. the rates for the semimajor axis 
and the eccentricity turn out to be zero; it is not so for the mean anomaly but no observational 
determinations exist for its extra-rate. 



2.1. Comparison with data 



According to the general outline of Section [H we may now consider a pair of planets A and 
B, take the ratio of their estimated extra-rates of perihelia and compare it to the prediction of eq. 
(|10p in order to see if they are equal within the errors 



Awt 



l g B (1 - e|) Zee 
a A (1 - el) \e A 



(11) 



If the modification of the gravitational potential not modelled by Pitjeva in estimating Am, 
accounts for what is unmodelled in the perihelia precessions, i.e. just for Am, then (jlip must be 
compatible with zero, within the errors. It must be noted that our approach is able to directly 
test the hypothesis that the proposed 1/r exotic force is not zero, irrespective of the magnitude of 
a because our prediction for the ratio of the extra-rates of perihelia, i.e. 

2 



Va 
Vb 



la B (l - e 2 A ) f e B 



-e B )\e A/ -1 + ^1- 
is just independent of a itself and is a specific function of the planets' orbital parameters a and e. 
Of course, the ratio of the perihelion rates cannot be used, by definition, to test the zero hypothesis 



■I + a/1 



(12) 



4 Indeed, V(r) of (47) in ( 



Adkins and McDonnell 



20071 ) is an additional potential energy. 
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Table 3: A B denotes the pair of planets used; II = Ari7 A /Ari7 B , A = (ob/oa) 1 / 2 an d 
F{exi ^b) is gi ven by eg. ffT^l) . The perihelion extra-rates for the inner planets have been 



retrieved from (IPitjeva 



2005bl ); their errors are not the 



extra-rates for the outer planets come from (IPitjeva 



orma 



, statistical ones. The perihelion 



2006al ); their formal errors have been 



re-scaled by a factor 10. T he uncertainties in the semimajor axes have been retrieved from 



Table IV of flPitieva 



2005a)): they are the formal ones, but, as can be noted, their impact is 
negligible. While II is always compatible with zero, this is definitely not the case for A. The 
eccentricity function F is always close to unity. 



A B 



n 



A 



F(eA,e B ) 



Mars Mercury 




0.03 ±0.2 





504 


± 


0(10" 




1 


008 


Mercury Jupiter 




0.6 ±4.1 


3 


666 


± 


0(10" 


- 9 ) 





989 


Earth Jupiter 




0.03 ±0.25 


2 


281 


± 


0(10- 


-10) 


1 


000 


Mars Jupiter 





02 ±0.17 


1 


847 


± 


0(10- 


-10) 





998 


Mercury Saturn 





004 ± 0.017 


4 


963 


± 


0(10- 


-9\ 





989 


Earth Saturn 





0002 ± 0.0011 


3 


088 


± 


0(10- 


-9\ 


1 


000 


Mars Saturn 




0.0001 ± 0.0009 


2 


501 


± 


0(10- 


-9 a 





998 


Jupiter Saturn 




0.006 ± 0.060 


1 


353 


± 


0(10- 


-9 a 


1 


000 


Mercury Uranus 




0.006 ±0.152 


7 


041 


± 


0(10- 


-8 a 





989 


Earth Uranus 




0.0003 ± 0.0087 


4 


380 


± 


0(10- 


-8 a 


1 


000 


Mars Uranus 





0002 ± 0.0048 


3 


459 


± 


0(10- 


-8a 





983 


Jupiter Uranus 





01 ±0.31 


1 


920 


± 


0(10- 


-8\ 





999 
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which, instead, can be checked by considering the apsidal precessions separately; indeed, in this 
case the prediction for the ratio of the perihelion precessions would be, by definition, 0/0. Note 
that, being Ave observational quantities, their ratio II = Acta/ Awe is a well defined quantity, 
with an associated uncertainty 511 which will be evaluated below. 
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From eq. (jlip and Table [3] it is possible to obtain 

^MarMer = 0.5 ± 0.2, 
*MerJup =4.2 ±4.1, 
^Earjup = 2.3 ± 0.2, 
^Marjup = 1.8 ± 0.2, 
^MerSat = 4.91 ± 0.02, 
*EarSat = 3.090 ± 0.001, 
^MarSat = 2.4984 ± 0.0008, 
^Jupsat = 1-36 ± 0.06, 

*MerUra = 6.9 ± 0.1, 
^EarUra = 4.383 ± 0.009, 
^MarUra = 3.543 ± 0.005, 
_ *J U pUra = 1-9 ± 0.3. 

The exact eccentricity-dependent factor 



F(e A ,e B ) 



<i_ e 2 ^ eB y /_i + A /i_ e 2 



1 - e| \e A J 



(13) 



-l + A /l-e| 



(14) 
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of eq. (jlip is always close to unity, so that its impact on the results of eq. (|13|) is negligible. The 
errors in ^ab due to 8a and 5 Am have been conservatively computed as 



Arir 



Aw B 



5Aw 



+ 



5Avo 



Aw A \ \Aw B 



+ 



lfa B Yf5a A 5a B \ 

+ 2 b U^ + ^ )F(eA;eB); 



(15) 



we did not optimistically summed the biased terms in a root-sum-square fashion because of 
the existing correlations among the estimated extra-precessions, although they are low with a 
maximum of aboul[] 20% between Mercury and the Earth. It turns out that by re-scaling the 



formal errors in the semimajor axes released by 



Pitieval (|2005al ) by a factor 10, 100 or more the 



results of eq. (|13p do not change because the major source of uncertainty is given by far by the 
errors in the perihelion precessions. In regards to them let us note that the tightest constraints 



Pitjeva 



(j2005al ) used only optical 



come from the pairs involving Saturn and Uranuqj for which 
data, contrary to JupiteJ^. Now, even by further re-scaling by 
perihelion extra-rates, i.e. by a factor 100 their formal, statistical uncertainties, the results do not 
change (apart from A= Jupiter B=Uranus). In the case of the pairs A=Earth, B=Uranus and 
A=Mars, B=Uranus even if the real error in the correction to the perihelion precession of Uranus 
was 1,000 times larger than the estimated formal one the answer of eq. (|13p would still remain 
negative. For A=Mars and B=Saturn a re-scaling of even a factor 10,000 of the formal error in 
the perihelion extra-rate of Saturn would be possible without changing the situation. In the case 
A=Earth, B=Saturn a re-scaling of 1,000 for the Saturn's perihelion extra-precession would not 
alter the result. 



5 L.I. thanks E.V. Pitjeva for such an information. 

6 But not for the pairs A=Uranus/Saturn, B=Saturn/Uranus yielding values of ^ab compatible 
with zero. 

7 This is why a re-scaling of 10 of the formal error in its perihelion extra-precession should be 
adequate. 
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As a conclusion, the logarithmic correction to the Newtonian potential of eq. ([!]) is ruled out 
by the present-day observational determinations of the Solar System's planetary motions. 



3. The power-law corrections in the f(R) extended theories of gravity 



In recent years there h as been a lot of interest in the so called f(R) theories of gravity 



(Sotiriou and Faraoni 



20081 ). In them, the gravitational Lagrangian depends on an arbitra ry 



analytic function / of the Ricci scalar curvature R (see (jCapozziello and Francaviglia 



20081 ) and 



references therein) . These theories are also referred to as "extended theories of gravity" , since 
they naturally generalize, on a geometric ground, GR, in the sense that when f(R) = R the action 



reduces to the usual Einstein- 



(jCapozziello and Francaviglia 



Rlbert one, and Einstein's theory is obtained. It has been showed 



2008 ) that these theories provide an alternative approach to solve, 
without the need of dark energy, some puzzles connected to the current cosmological observations 
and, furthermore, they can explain the dyn amics of the rotat i on cu rves of the galaxies, without 



requiring dark matter. Indeed, for instance 



Capozziello et al 



d2QQ2l), starting from f(R) = f R k 



obtained a power-law correction to the Newtonian gravitational potential of the form 



P 



GM ( r 



(16) 



(where (5 is related to the exponent k of the Ricci scalar R), and they applied eq. (fT6|) to a 
sample of 15 low surface brightness (LSB) galaxies with combined HI and Ha measurements of 
the rotation curve extending in the putative dark matter dominated region. They obtained a very 
good agreement between the theoretical rotation curves and the data using only stellar disk and 
interstellar gas when the slope k of the gravity Lagrangian is set to the value k = 3.5 (giving 
f3 = 0.817) obtained by fitting the SNela Hubble diagram with the assumed power-law f(R) model 
and no dark matter. 



Here we wish to put on the test eq. (|16p in the Solar System with the approach previously 
examined. 
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First, let us note that an extra-radial acceleration 

(J3 - 1)GM 



(17) 



can be obtained from eq. f)16[> ; let us now work out the secular precession of the pericentre of a 
test particle induced by eq. (fT7j) in the case (3 — 2 < 0. By proceeding as in Section [2] we gee 

(P-1)VGM _3 



2vrr£ 



■■■G(e-(3), 



with 



2- 



cos E — e 
(1-ecos^) 2 "/ 3 



dE. 



(18) 



(19) 



Since we are interested in taking the ratios of the perihelia, there is no need to exactly compute eq. 
(I19|) : the eccentricities of the Solar System planets are small and similar for all of them, so that 
we will reasonably assume that G{ex-n)/G{e-Q] (3) ~ 1. Note that|j (3 = would yield a vanishing 
apsidal precession because G(e; 0) = 0. In fact, the case j3 = is compatible with all the estimated 



extra-rates of Table [T] and Table [2j (3 = 0, within t 



performed in the Solar System by 



re errors, is also the outcome of different tests 



Zakharov et al 



(|2006h . 



3.1. Comparison with data 



According to 



Capozziello et al 



(|2007l ). a value of j3 = 0.817 in eq. (fT6|) gives very good 
agreement between the theoretical rotation curves and the data, without need of dark matter. 
Again, if we consider a pair of planets A and B and take the ratio of their estimated extra-rates of 



perihelia, by taking f3 = 0.817, we may define 

Acta 



OB 
OA 



0.683 



(20) 



Adkins and McDonnel] 



20071 ) 



See also 

9 The case (3 = 1 is not relevant because it would yield a constant additional potential and no 
extra-force. 
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Of course, such kind of test could not be applied to the (3 = case since, in this case, we 
would have the meaningless prediction Va/Vb = 0/0. If the modification of the gravitational 
potential of eq. (|16f) exists and is accounted for by the estimated corrections Azu to the standard 
Newton-Einstein perihelion precessions, then the quantities Hab must be compatible with zero, 
within the errors. Instead, what we obtain, from Table [U Table [2j Tableland eq. (|18|) . is 
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Table 4: A B denotes the pair of planets used; II = Ari7 A / 'Awb, B = 
(ob/oa) ' 683 (corresponding to = 0.817). The perihelion extra-rates for the inner plan- 



ets have been retrieved from (jPitjeva 



2005bl ); their errors are not the forma l, statistical ones. 



The perihelion extra-rates for the outer planets come from (jPitjeva 



2006al ); their formal er- 



rors have been re-scaled by a factor 10. T he uncertainties in the semimajor axes have been 



retrieved from Table IV of (jPitjeva 



2005a ): they are the formal ones, but, as can be noted, 
their impact is negligible. While II is always compatible with zero, this is definitely not the 
case for B. 



A B 


n 


B 






Earth Mercury 


0.05± 0.19 


0.523 


± 


0(10- 


Mars Mercury 


-0.03 ±0.2 


0.392 


± 


0(10" 


Mercury Jupiter 


-0.6 ±4.1 


5.898 


± 


0(10" 


Earth Jupiter 


-0.03 ±0.25 


3.084 


± 


0(10" 


Mars Jupiter 


0.02 ±0.17 


2.313 


± 


0(10- 


Mercury Saturn 


0.004 ± 0.017 


8.921 


± 


0(10- 


Earth Saturn 


0.0002 ± 0.0011 


4.665 


± 


0(10- 


Mars Saturn 


-0.0001 ± 0.0009 


3.499 


± 


0(10- 


Jupiter Saturn 


-0.006 ± 0.060 


1.512 


± 


0(10- 


Mercury Uranus 


-0.006 ±0.152 


14.383 ±O(10 


Earth Uranus 


-0.0003 ± 0.0087 


7.522 


± 


0(10- 


Mars Uranus 


0.0002 ± 0.0048 


5.642 


± 


0(10- 


Jupiter Uranus 


0.01 ±0.31 


2.438 


± 


0(10- 



13 ) 

13^ 
9 
9 
10 
8 
9 
9 
9 
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SEarMer = 0.5 ± 0.2, 
^MarMer = 0.4 ± 0.2, 
^MerJup = 6.5 ± 4.2, 
^EarJup = 3.1 ± 0.2, 
^MarJup = 2.3 ± 0.2, 
^MerSat = 8.92 ± 0.02, 
< ^EarSat = 4.666 ±0.001, (21) 
^MarSat = 3.4997 ± 0.0008, 
HjupSat = 1-52 ± 0.06, 

S M crUra= 14.4 ±0.1, 
S Ea rUra = 7.523 ± 0.008, 
^MarUra = 5.642 ± 0.005, 
^JupUra = 2.4 ± 0.3. 

It is remarkable to note that even if the formal error in the Saturn's apsidal extra-precession was 
re-scaled by a factor 10,000 instead of 10, as done in this paper, the pairs A=Earth B=Saturn and 



-20- 



A=Mars B=Saturn would still rule out eq. (|18p . A re-scaling by 1,000 of the Uranus' perihelion 
extra-precession would still be fatal, as shown by the pairs A=Earth B=Uranus and A=Mars 
B= Uranus. 



Thus, also the power-law correction to the Newtonian potential of eq. (|16|) with f3 = 0.817 
is ruled out. Criticisms to R k models of modified gravity were raised on different grounds 



by 



Noiiri and Qdintsovl (|2003l ) who proposed more realistic models in (|Cognola et al 



2008 ) 



Moreover, it is well-established now that DM shows also particle-like properties. In this respect, 
the proposal of R k gravity as DM (thanks to a change of the Newton potential) is not considered 
as a realistic one now. A more realistic DM candidate from f(R) gravity was suggested by 



Nojiri and Qdintsovl (|2008l ). They show that not only a correction to the Newton potential 
appears, but also composite graviton degree of freedom shows particle-like behavior, as requested 
by DM data. 



4. Comments and conclusions 

In this paper we have studied the secular precession of the pericentre of a test particle in 
motion around a central mass M whose Newtonian gravitational potential exhibits a correction 
which has a logarithmic and power-law behavior. In order to put on the test the hypothesis that 
such extra-forces are not zero we devised a suitable test by taking into account the ratios (and 
not the extra-rates of the perihelia of each planet at a time separately, or a linear combination 
of them, since their uncertainties would fatally prevent to obtain any useful constraints) of the 
corrections to the secular precessions of the longitudes of perihelia estimated by E.V. Pitjeva for 
several pairs of planets in the Solar System. The results obtained, resumed by eq. (|13|) and eq. 
(|2ip . show that modifications of the Newtonian potentials like those examined in this paper are 
not compatible with the currently available apsidal extra-precessions of the Solar System planets. 
Moreover, the hypothesis that the examined exotic force terms are zero, which cannot be tested 
by definition with our approach, is compatible with each perihelion extra-rate separately, in 
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agreement with our results. It must be noted that to give the hypothesized modifications of the 
Newtonian law the benefit of the doubt, and given that the formal errors in a and Am of the outer 
planets are probably underestimates of the true uncertainties, we multiplied them by a factor 
ten or even more, as suggested to one of us (L.I.) by some leading experts in the ephemerides 
generation field like E.M. Standish, but the answers we obtained were still negative. Another 
thing that should be pointed out is that, in principle, in assessing 'I'ab and Hab one should have 
used Am* = Am — ^canonical; where (^canonical represents the mismodelled part of the modelled 
standard Newton/Einstein precessions. However, neglecting them did not affect our conclusions, 
as can be easily noted by looking at the II and A/B columns in Tableland Table HI Indeed, the 
residual preces sions due to the imperfect knowledge of, e.g., th e solar quad rupole mass moment 



(jPitieva 



2005b] ) 5.J2/J2 ~ 10% are of the order of < 10~ 3 " cy- Mlorio 



2007d ). and even smaller are 



the mismodelled pre cessions due to other potential sources of errors like the asteroid ring or the 



Kuiper Belt objects (jlorio 



2007cf ) 



However, caution is in order because, at present, no other teams of astronomers have 
estimated their own corrections to the Newtonian/Einsteinian planetary perihelion rates, as it 
would be highly desirable. If and when we will have, say, two independent determinations of the 
anomalous perihelion rate of a given planet x = x^ est ± 5x and y = Ubest ± Sy we could see if they 
are compatible with each other and take their difference latest — 2/best| as representative of the real 
uncertainty affecting the apsidal extra-precession of that planet. Moreover, it would be interesting 
to see if for different sets of estimated corrections to the perihelion rates the figures for n in Table 
[3] and Table H] change by an extent sufficient to alter the conclusions of eq. f|13j) and eq. (|2ip . 
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